Massively parallel hybrid algorithm for simulation of seismic wave propagation in complex media

This paper presents a problem-oriented approach, designed for numerical simulation of seismic wave propagation in models containing geological formations with complex properties such as anisotropy, attenuation, and small-scale heterogeneities. Each of the named property requires special treatment which increases computational complexity of the algorithm in comparison with ideally elastic isotropic media. At the same time these formations are typically relatively small taking up to 25% of the model, thus the computationally expensive approaches can be used only locally to speed-up the simulations. In this paper we discuss both mathematical and numerical aspects of hybrid algorithm paying the main attention to its parallel implementation.